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ABSTRACT 


Stress  propagation  in  a  nonlinear  viscoelastic  material 
is  analyzed  for  the  unloading  portion  of  a  complete  stress- 
strain  cycle.  The  loading  process  is  the  subject  of  a 
previous  report.  The  elastic  unloading  theory  developed 
by  Malvern  is  reviewed  and  extended  to  include  stress-strain 
laws  with  nonlinear  unloading.  Numerical  techniques  for 
solution  are  outlined  for  both  infinitesimal  and  finite 


strain  theories 


LIST  OF  SYMBOLS 


cr  - 

V  - 
£  - 

oi  - 

E  - 

at-j\r  - 

f  - 
c  . 

±  . 

9C  - 


stress 

particle  velocity 

infinitesimal  strain  or  space  gradient  of  displacement 
finite  strain  = 

Young ' s  modulus 

dynamic  constants  of  material 

density 

velocity  of  initial  wave  front 
time 

Lagrangian  space  variable 


All  other  terms  are  defined  in  the  text  where  first  introduced. 


INTRODUCTION 


Experimental  evidence  from  the  dynamic  testing  of  metals  and 

eafth  materials  indicates  that  upon  removal  of  external  forces 

exceeding  the  yield  point,  the  material  does  not  return  to  its 

initial  state.  A  commonly  accepted  procedure  for  dealing  with 

this  hysteresis  effect  is  to  assume  that  the  stress-strain  laws 

in  the  loading  and  unloading  cycles  are  not  the  same.  The  material 

is  often  assumed  to  unload  elastically  with  slope  equal  to  that 

of  the  initial  elastic  loading  path.  Boundary  value  problems 

associated  with  materials  that  behave  plastically  on  loading  and 

1 

elastically  on  unloading  have  been  discussed  by  Lee.  Materials 

that  behave  plastically  on  both  loading  and  unloading  have  been 

2 

investigated  by  Norland.  Unloading  phenomena  in  strain-rate 

.  .  .  3 

sensitive  materials  have  been  investagated  by  Malvern  with  a 
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subsequent  contribution  by  Rubin.  Both  Malvern  and  Rubin  con¬ 
fine  their  attention  to  elastic  unloading  parallel  to  the  initial 
slope  of  the  static  stress-strain  curve. 

The  mathematical  analysis  is  simplest  when  the  characteristic 
directions  are  unchanged  during  the  transition  from  loading  to 
unloading.  This  suggests  the  use  of  infinitesimal  strain  theory 
and  a  linear  unloading  law  such  that  unloading  shock  waves  propa¬ 
gate  with  constant  velocity  equal  to  that  of  the  initial  wave 
front.  Under  these  conditions  the  mechanical  shock  conditions 
are  sufficient  to  determine  the  subsequent  motion.  In  the  general 


case  firhen  the  characteristic  directions  in  front  of  and  behind  the 
shock  are  not  parallel,  we  must  assume  some  form  for  the  dependence 
of  internal  energy  on  stress  and  strain  which  is  utilized  in  form¬ 
ing  the  Hugbniot  function.  The  Hugoniot  function,  in  conjunction 
with  the  mechanical  shock  conditions,  is  sufficient  to  determine 
the  jumps  in  stress  and  strain  across  the  unloading  shock.  Since 
the  stress-strain  law  that  the  material  obeys  once  unloading 
commences  differs  from  that  during  loading,  it  is  plausible  also 
to  use  a  different  internal  energy  function.  However,  we  will 
show  that  this  assumption  prevents  the  shock  wave  from  ever  being 
cotnpletely  absorbed.  The  situation  for  different  internal  energy 
functions  across  the  shock  is  similar  to  that  for  detonation  waves 
where  an  explosive  mixture  undergoes  a  change  of  state  across  the 
detonation  front.  Processes  of  this  type  are  discussed  by  Courant 
and  Friedrichs^;  however,  their  analysis  cannot  be  directly  applied 

to  the  problem  at  hand  unless  we  generalize  the  concept  of  the 
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Hugoniot  function  along  the  lines  suggested  by  Duvall.  This 
generalization  takes  into  account  the  region  of  variable  state 
into  which  the  shock  wave  moves. 

Because  of  the  intractability  of  the  mathematics  in  the  general 
case,  the  physical  processes  are  obscured.  To  obtain  some  insight, 
into  the  problem,  we  review  Malvern ' s  description  of  the  elastic 
unloading  process  from  a  more  descriptive  point  of  view  before 
considering  the  general  case. 
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Malvern's  discussion  of  the  unloading  phenomena  when  the 
stress  at  X  s  0  is  instantaneously  reduced  to  zero  and  main¬ 
tained  at  that,  value  follows.  See  Figure  I.  An  unloading 
shock  wave  propagates  along  the  characteristic  from  the 
point  X  s  0,  t  =  t^  at  Srhich  the  stress  is  released.  This 
shock  decreases  in  amplitude  and  may  or  may  not  be  absorbed. 
The  boundary  between  the  elastic  and  viscoelastic  states  is 
identical  with  the  shock  path  until  the  shock  strength  goes 
to  zero,  at  which  point  the  elastic-viscoelastic  boundary 
leaves  the  characteristic  and  must  be  found  by  a  trial- 
and-error  method.  Rubin,  after  linearizing  the  equations, 
uses  asymptotic  methods  to  show  that  an  unloading  shock  is 
always  absorbed  in  a  finite  time. 

We  present  an  alternative  method  for  calculating  the 
elastic-viscoelastic  boundary  after  it  leaves  the  charac¬ 
teristic  and  its  subsequent  interaction  with  elastic  waves. 
While  this  alternative  method  appears  to  be  more  adaptable 
to  programming  for  numerical  calculation  than  Malvern's, 
the  two  methods  are  equivalent  and  should  give  identical 
answers.  Malvern's  description  of  the  unloading  phenomena 
is  also  modified  to  include  the  case  that  the  material  behind 
the  unloading  shock  is  in  the  viscoelastic  state.  This  theory 


is  illustrated  by  tracing  qualitatively  the  changes  a 

rectangular  stress  pulse  of  amplitude  01  and  duration  t 

o  O 

applied  at  x  =  0  undergoes  as  it  propagates  into  the  material. 


THEORY 


The  dynamic  behavior  of  the  material  can  be  illustrated  by 
outlining  those  regions  of  the  stress-strain  plane  that  contain 
physically  attainable  values  for  these  variables.  In  Figure  II, 

9. CO")  ^  is  the  static  stress-strain  law  during  loading. 

The  line  0~  —  B  S.  is  the  continuation  of  the  initial  slope 

of  the  curve  g_((r):=  -  The  line  0~  ~  ^  is  the 

locus  of  stress-strain  pairs  at  infinite  strain  rate.  Regions 
II  and  III  are  the  locus  of  points  for  the  viscoelastic  state. 
Region  I  is  the  locus  of  points  where  the  material  has  become 
elastic;  any  point  in  Region  I  must  be  connected  to  a  point  in 
Regions  II  or  III  by  an  unloading  shock  wave. 

A  material  governed  by  Malvern's  theory  could  only  assume 
values  in  Region  II  of  Figure  II  during  the  loading  portion  of 
a  complete  stress  cycle.  Hence,  Malvern's  theory  predicts  that 
Young's  modulus  E  calculated  from  any  dynamic  test  will  never 
exceed  E  obtained  during  a  static  test.  However,  experimenters 
in  soil  mechanics  report  increases  as  high  as  30!^  in  Young's 
modulus  when  measured  dynamically  over  the  value  obtained  in 
static  tests. ^  The  present  theory  adds  Region  III  to  the  set 
of  attainable  stress-strain  pairs  of  Malvern's  theory  to  allow 
this  apparent  increase  in  Young's  modulus  with  increasing 
strain  rate.  This  generalization  of  Malvern's  theory  is 
obtained  by  introducing  an  additional  dynamic  constant  b  to 


supplement  Malvern's  dynamic  constant  k  which  is  identical  to 
l/a  in  the  present  theory. 

Qualitatively  correct  wave  shapes  can  be  deduced  without 
resorting  to  numerical  analysis  for  the  special  case  of  in¬ 
finitesimal  strain  theory  and  linear  unloading  with  slope  b/a 
from  any  point  in  Region  I.  Under  these  assumptions,  the 
characteristic  directions  for  the  complete  stress  cycle  are 
the  family  of  straight  lines  X—  t  f +  constant. 

In  contrast  to  plastic  materials,  stress  and  strain  are  not 
constant  along  these  characteristics  because  nonlinear  visco¬ 
elastic  materials  do  not  possess  simple  wave  solutions.  A 
physical  explanation  for  the  absence  of  a  simple  wave  solu¬ 
tion  is  discussed  in  a  previous  report.®  However,  much 
physical  insight  can  still  be  gleaned  by  considering  the  vari¬ 
ation  of  stress  and  strain  along  individual  character¬ 

istics  as  a  function  of  time.  With  reference  to  Figure  III  , 
<U”  and  ^  along  are  related  by  U~  ~  ^  E. 

This  is  the  line  segment  A-A’  in  the  stress-strain  plane 
shown  in  Figure  IV.  The  origin  in  Figure  III  where  the  stress 
rises  discontinuously  from  zero  to  transforms  to  the 

point  A  in  Figure  IV.  As  this  decaying  shock  propagates  into 
the  material,  the  associated  values  of  0“  and  £_  move 
down  the  line  segment  from  A  to  A'  as  t,  increases  from 

O  to  oo 
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The  values  to  stress  and  strain  associated  with  the  point 
B,  which  is  where  ^  intersects  x=  O,  differ  from  those 
associated  with  A  because  the  strain  has  increased  in  the  time 


interval  from  A  to  B.  The  actual  path  taken  through  the  0~j  £, 
plane  along  is  not  known  but,  because  of  continuity,  it 

must  be  reasonably  close  to  A-A' ,  never  intersecting  A-A' o 
This  path  is  shown  schematically  as  B-B' .  Thus  the  motion 

can  be  described  qualitatively  by  plotting  the  variation  of 

( ■'1 ) 

(J~  and  along  any  C^.  characteristic  through 

the  £  plane  from  at  x  =  0  to  some  point  on  ^ 

as  t  approaches  OO  .  These  paths  through  the  CTj  £ 

plane  are  shown  as  A-A',  B-B',  C-C,  D-D',  ecto,  in  Figure  IV. 

(n)  .  .  .  . 

At  t^,  where  C+  originates,  the  stress  is  instantan¬ 
eously  reduced  to  zero.  Then  the  values  of  stress  and  strain 
at  the  back  of  the  resulting  shock  are  connected  to  the  values 
in  front  of  the  shock,  which  are  plotted  at  the  point  ^ 
by  a  line  parellel  to  A-A',  terminating  at  CT  —  O  •  i.e,, 
the  point  ^2  in  Figure  IV.  The  values  of  stress  and 
strain  in  front  of  and  in  back  of  the  shock  are  connected  by 
this  path  because  unloading  in  the  viscoelastic  Regions  II 
and  III  follows  such  paths  for  the  infinite  strain  rate  gen¬ 
erated  by  the  step  in  stress.  By  hypothesis,  the  material 
unloads  in  Region  I  along  this  path,  independently  of  the 


strain  rate. 
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At  X  8  0,  we  were  able  to  calculate  £.  across  the 
shock  because  we  knew  A  0~  from  the  boundary  condition. 

However  when  x  0,  we  will  have  to  calculate  ^  0~  from  the 
characteristic  equations  and  the  mechanical  jump  conditions. 
Note  that  the  foot  of  the  shock  is  well  within  Region  I,  at 
least  for  small  x.  This  situation  is  identical  to  Malvern's, 
and  the  equation  for  A  ^  in  tefms  of  the  known  quantities 
along  O,  -  is  ZrT.o 

ACr(§,t-t,J=  Acrfo,o)-5^  (1) 

^  o 

This  equation  is  derived*  by  subtracting  the  characteristic 
equations  valid  along  in  front  of  and  behind  the  shock 

and  using  the  mechanical  jump  conditions. 

Rubin  showed  that  a  linearized  form  of  (1)  always  van¬ 
ishes  for  some  finite  value  of  t,  say  point  Q  in  Figures  III 
and  IV,  Before  (1)  vanishes,  will  reach  a  point  P 

such  that  the  values  of  stress  and  strain  at  the  foot  of  the 
shock  lie  on  ^  "fC^)  >  i-e.,  the 

point  in  Figure  IV.  Beyond  this  point  we  assume  that 

the  material  is  viscoelastic  since  the  foot  of  the  shock  lies 
above  the  static  curve.  Because  of  this  change,  (1)  is  no 
longer  applicable  for  computing  ^  0~  .  With  the  mechanical 

jump  conditions  and  the  altered  characteristic  equations 
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*See  Appendix  II. 


in  front  of  and  behind  the  shock,  (1)  becomes* 


r 


where 


ai^£,  = 

stress , 

strain  along 

P-  a-  a, 

£j.= 

= 

stress , 

time  at 

strain  along 

point  P 

AcrCP”)  = 

ACT 

at  P  calculated  from 

The  stress-strain  path  at  the  foot  of  the  shock  is  no 
longer  ^  shown,  but  follows  some  path 

<4-  (4'  -  <i,  with  a  discontinuity  in  slope  at  0^. 

This  can -be  seen  by  rewriting  (1)  and  (2)  in  differential 
form  and  sketching  their  graphs.  Equations  (1)  and  (2) 
can  be  written  as 


o/idCT 

J  t 


[  ^((TT)  - -f  (•£,)  J 


(3) 


ci^  _  _/ 
Jtr  ^ 


( po:)-^(oi))  -  -  fa^)J  (4) 


where  subscripts  1  and  2  refer  to  points  along  the  paths 
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*See  Appendix  II. 


and  ^  •  In  (3)  the  term  -  fCi.,)  , 

at  a  typical  point  t<tp  ,  is  represented  graphically  as  the 

line  segment  R-R  .  At  time  t^,  •^(^)“  'f(£i') 

represented  graphically  as  the  line  segment  —  oC,  ; 

from  Figure  IV,  it  is  a  monotonically  decreasing  function  of 

time  as  R  moves  from  to  If  the  material  did  not 

cross  the  elastic-viscoelastic  boundary,  R  would  continue 

along  0i  —  P  to  Q  where  the  shock  would  be  absorbed 

and  then  from  Q  to  where  •*  ^(^i) 

would  go  to  zero  at  t  =  oO  .  Hence  we  can  visualize  (3) 

graphically  as  shown  in  Figure  V. 

From  (1)  we  see  that  the  area  under  this  curve  from 
X  a  0  to  tg  is  numerically  equal  to  (OjO) 

Equations  (3)  and  (4)  are  equal  t  ,  because  at  this 

point  the  second  term  on  the  right  of  (4)  is  zero.  From 
Figure  IV  it  is  apparent  that 

for  all  t  and  that  pai)--f(£,)  ^  f  ( £  ^) 

for  t  t  p  where  the  foot  of  the  shock  lies  in  Region  II. 

Thus  for  t  >  t  p  ,  computed  from  (4)  will  be 


’  JT 


computed 


below  the  corresponding  value  computed  from  (3).  We  can 
also  show  that  at  t  =  t  „  , 

P 

from  (3)  and  (4)  differs  by  an  amount 

ciw  \  4!^']  ^ 

=  -ZCL. 
tp 


Jt’-  (H) 
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Therefore,  in  Figure  V  there  will  be  a  discontinuity  in  slope 
at  t  p  ;  the  new  locus  of  points  will  plot  below  the  curve 
computed  from  (3)  for  all  t  t  p  as  shown  in  Figure  VI.  The 
time  at  which  the  shock  is  absorbed  must  be  increased  to  tg^^  ^ 
so  that  the  area  under  the  new  curves  in  Figure  VI  from  t  p 
to  tg^  is  equal  to  the  area  under  the  curve  from  t  p  to  t^ 
in  Figure  V.  We  must  now  admit  the  possibility  that  the  shock 
is  never  absorbed.  From  Figure  IV  ,  the  stress-strain  path 
at  the  foot  of  the  shock  now  branches  at  with  a 

discontinuity  in  slope,  and,  if  absorbed,  it  intersects  the 
path  in  front  of  the  shock  at  the  point  The  path  seg¬ 
ment  ^  is  shown  intersecting  the  path 

with  different  slope.  Although  this  is  not  proved  and  must 
be  verified  by  numerical  calculations,  it  seems  physically 
plausible  that  when  the  shock  strength  goes  to  zero,  a  dis¬ 
continuity  in  the  derivatives  of  stress  and  strain  remains. 
This  can  be  visualized  by  realizing  that  the  peaks  of  the 
resulting  wave  forms  are  cusped.  If  this  discontinuity  did 
not  exist,  the  wave  peaks  would  be  smooth. 

As  shown  in  Figure  III,  the  x-t  plane  can  be  divided  into 
four  areas.  Area  E  is  the  loading  pulse  and  is  characterized 
by  increasing  stress  and  strain  as  functions  of  tiroe«  These 
stress-  strain  points  fall  in  Regions  II  and  III  of  Figure  II; 
the  stress-strain  paths  are  the  lines  A-A' ,  B-B' ,  etc^,  in 
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I 

I 

1 

1 


Figure  IV.  Area  F  is  characterized  by  being  entirely  elastic 
with  the  stress,  strain,  and  particle  veldcity  known  along 
the  foot  of  the  unloading  shock  which  precedes  it.  Area  G 
is  also  elastic  but  differs  from  Area  F  because  the  visco¬ 
elastic-elastic  boundary  preceding  it  is  unknown  and  must 
be  calculated  by  special  techniques.  Although  Area  H  is 
viscoelastic,  it  differs  from  Area  E  because  the  stress 
and  strain  are  decreasing  functions  of  time. 

The  solution  in  Area  F  can  be  obtained  by  the  principle 
of  superposition.  Assume  that  the  necessary  calculations 
have  been  made  to  obtain  ,  E.  ,  and  [T  at  the  foot  of 
the  unloading  shock.  Let  S  be  a  parameter  measuring  dis¬ 
tance  along  the  unloading  shock.  Along  this  path  Q”  , 

,  and  \r  are  given  by 

(jcs)  ^  (j;  " 


£(S)  - 

V  (  S )  r  \/,  f  ^ 

Along  the  boundary  x  =  0  the  stress  is  O  for  t  ^ 
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By  tbs  principle  tif  superposition,  the  stress  st  any  point 

(K,t)  in  Area  P  is 

(TfX.t)  =  f-  i  0-(i)Js  -  J  O-^rs)  0/  S  <=) 

%  0 

The  stressx'strain  relation^  in  Area  F  is 

crrsj  -  0"  =  (6) 


with  i'5'),  (6)  becoines 


r 


£(x,tJ=  f  t  0~U)J5-t  CTfsJ  J5 


Note  that  =  E  in  Area  F. 

We  now  wish  to  trace  the  elastic-viscoelastic  boun¬ 
dary  after  it  leaves  the  characteristic »  Moving 

along  the  characteristic  from  x  »  O,  we  first 

traverse  the  elastic  Area  F  where  the  solution  is  known. 
This  is  the  path  segment  in  Figure  IV 

where  (f?^  represents  the  boundary  between  Areas  F 

and  G. 


The  calculation  of  the  elastic-viscoelastic  boun¬ 
dary  will  be  illustrated  in  Figure  VII  which  is  an  en¬ 
largement  of  Figure  HI  in  the  neighborhood  of  P.  We 
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know  that  -  '^(^)  along  thf  elattlc-viscoelastic 

boundary  after  it  leaves  the  characteristic.  At 

point  Where  this  boundary  cuts  ,  the  follow¬ 

ing  characteristic  equations  hold: 

a.  Along  C_ 

b.  Along 

(o-p,  jrO-f  c"  -fca.,)] /iz 


At  CT  ,  £  .  ,  and  ]/~  are  related  by  the  charac¬ 

teristic  equations  for  a  non-linear  elastic  material  governed 


by  ^(cr)  --((£) 

£p 

V,  -f 

where  ^  and  C  C  O' ) 

fiz)  =  ^(^) 

(9)  can  be  written 


I 

^  Jcr 


J 


C(^J 


(10) 


are  derivable  from 

With  (10),  (8)  and 


(C^-0^^)4-^C  -  2.0^  -f(’£ciO(ll) 

L  J  o  J  ^ 
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(12) 


( <^p,  -  '^=.'1  -  ih<T)  ■  S " '^■'> 

which  gives  0^  implicitly  as  function  of  known  quantities 
(n) 

along  .  The  calculation  procedure  follows:  aseume 

'X’P^  =  'Xp  s  solve  (11)  and  (12)  for  (j^  ,  and  compare 

them.  If  the  values  of  Q”p  calculated  from  (11)  and  (12) 

-V/  (n+l)  ^  c 

are  unequal,  move  /C.'yj  along  an  increment  J 

and  repeat  the  calculation.  Repeat  this  procedure  until 

the  S'  from  (11)  and  (12)  are  equal.  This 

is  the  correct  location  of  the  boundary. 

This  procedure  for  calculating  the  elastic-viscoelastic 

boundary  assumes  no  interaction  with  other  waves.  However, 

as  soon  as  the  viscoelastic  boundary  leaves  the  C  charac- 

teristic,  elastic  unloading  waves  overtake  it  and  modify 

both  its  position  and  associated  values  of  stress,  strain, 

and  particle  velocity. 

To  simplify  the  analysis,  consider  the  unloading 
wave  as  made  up  of  TL  stress  increments  of  magnitude 
^  Cr  separated  by  time  intervals  A  T.  .  The 

first  unloading  increment  propagates  along 
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through  th«  clastic  Region  P.  There  is  no  interaction  in 
this  region  except  a  change  in  the  level  of  the  stress  incre¬ 
ment  as  it  encounters  the  rising  stress  level  radiated  from 
the  back  of  the  unloading  shock.  This  unloading  increment 
meets  the  elastic-viscoelastic  boundary  at  the  point  in 
Figure  VII.  Since  the  strain-rate  associated  with  this 
step  is  infinite,  changes  in  stress  and  strain  are  given  by 
the  shock  conditions.  Hence  at  we  know  both  AU' 
and  A  ■>  This  situation  is  now  identical  to  the  previous 

( n )  .  .  I  — 

analysis  along  C  at  the  origin.  The  increment 

(n+l)  .  . 

will  move  into  the  viscoelastic  area  along  .  Initially, 

the  foot  of  this  small  shock  will  be  below  the  elastic- 
viscoelastic  boundary.  The  decrease  in  magnitude  of  the 
increment  A  Q~  ,  as  it  penetrates  the  viscoelastic 
area,  is  governed  by  (1)  until  a  new  point  P  along 
is  reached,  such  that  the  foot  of  this  shock  lies  on 
^(0-)  ^  -fra)  .  At  this  point  the  elastic-viscoelastic 
boundary  again  leaves  the  characteristic;  the  law  con¬ 

trolling  the  decay  in  shock  strength  changes  from  (1)  to 
(2).  Further  along  ,  the  increment  may  eventually 

be  absorbed  at  a  point  Q^.  The  n  unloading  stress  increments 
impinge  upon  the  elAstic-visceelfiStio  boundary  in  succession 
along  the  Shaxacteristic  as  shown  in  Figure  VIII.  The 
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procedure.-can  be  extended  indefinitely  by  generalizing  the 
parameter  S  and  hence  the  functions  Q~C^)  ,  etc.,  in  (5) 
and  (7), so  that  they  follow  the  elastic-viscoelastic  boundary 
after  it  leaves  the  characteristic  at  P.  For  this 

extension  the  limits  of  integration  must  be  modified.  This 
procedure  extends  the  region  in  which  (5)  and  (7)  are  valid 
to  include  G  in  Figure  III. 

The  curvature  of  the  boundary  after  it  leaves  can 

also  be  deduced  from  this  analysis.  As  shown  in  Figures  III, 
VII,  and  VIII,  the  velocity  associated  with  the  boundary 
decreases  with  time  and  distance.  This  immediately  raises 
the  question  whether  the  boundary  even  comes  to  rest,  so 
that  some  portion  of  the  material  remains  viscoelastic 
and  hence  deforms  forever.  This  will  be  shown  not  to  be 
the  case.  The  velocity  of  the  boundary  approaches  a  limit 
controlled  by  the  static  parameters  of  the  material.  Re¬ 
call  that  along  this  boundary  C(r)  =  and  that 

changes  in  stress,  strain,  and  particle  velocity  are  inter¬ 
related  by  equations  (10)  which  are  derived  from  the  charac¬ 
teristic  equations  of  a  plastic  material.  Hence  along 
this  boundary  -  and  only  this  boundary  -  the  material  behaves 
plastically  with  the  stress-strain  curve  -p 
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Thess  characteristic  equations  are  known  to  be  valid  only 
along  the  characteristic  directions  associated  with  the 
stress  and  strain  at  a  given  point.  Since  and  £.  vary 
along  this  boundary,  the  characteristic  direction  con¬ 
tinuously  varies;  hence  the  boundary  must  be  an  envelope 
of  characteristic  directions  of  the  plastic  material 
defined  by  the  static  stress-strain  curve.  Since  the  stress 
and  strain  decrease  along  this  boundary,  the  associated 
velocity  will  increase  or  decrease  as  the  static  curve  is 
concave  or  convex  to  the  strain  axis.  The  limiting 
velocity  which  the  elastic-viscoelastic  boundary  approaches 
is  determined  by  the  initial  slope  of  ^  Ccr  )  =  'f  (L)  ; 

i.e.,  C  5=  ^  ^  ^  ^  »  which  must  be  less  than  or  equal  to 

f  Ir- 

C  =  VcC^  •  Regardless  of  whether  the  velocity 

increases  or  decreases,  the  material  becomes  entirely 

T  .  rt 

elastic.  However  if  |  W  ,  the  viscoelastic 

region  increases  in  width;  while  if  they  are  equal,  the 
width  approaches  a  limiting  value. 

Examination  of  (5)  alnd  (7)  shows  that  as  t— ?  00  , 

the  stress  distribution  goes  to  zero  leaving  a  non-zero 
strain  distribution  calculated  from 


£(><^00)=  £(S)  -  ^  CTfs) 


(13) 
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We  are  now  able  to  describe  qualitatively  the  result¬ 
ing  Wave  shapes,  having  considered  how  CT*  and  S. 
vary  along  the  C^.  family  of  characteristics  and  noting  the 
changes  required  by  conservation  of  mass  and  momentum  as 
the  boundaries  between  various  regions  in  the  x-t  plane 
are  crossed. 

First,  consider  the  time  variation  of  0~  and  £1 
at  X  =  0.  The  former  is  known  since  it  is  the  boundary 
condition.  The  latter  is  known  to  be  an  immediate  elastic 
response,  ,  followed  by  flow  asymptotically  tending 

to  the  value  .  The  quantity  E.  eso  given  by 

the  solution  of  ^  (  0^  ]  r  C  ^  ooj 

is  shown  in  Figure  II.  Hence  Eoo  is  the  maximum 

strain  .the  material  would  attain  if  there  were  no  time- 


rate  effects, with  the  same  boundary  condition.  At  t^, 
the  stress  is  instantaneously  reduced  to  zero.  This  is 
illustrated  by  the  path  (p2^  in  Figure  IV. 

A  sketch  of  these  wave  shapes  is  shown  in  Figure  IX. 

Consider  the  time  variation  of  ^  and  £.  at 


a  point  intermediate  between  x  =  0  and  x  p  in  Figi 

The  stress  remains  zero  until  a  time  t  =  x/c  when  the 


in  Figure  III. 


initially  discontinuous  wave  fronts  of  amplitude 


crC4)  . 


arrive.  These  wave  fronts  are  followed  by  a 
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I 

I 

I 


period  of  increasing  stress  and  strain  tending  asynptot- 
ically  to  the  values  until  a 

time  't-o  f  '^/c  when  the  unloading  shock  wave 

arrives , instantaneously  reducing  the  stress  and  strain. 
The  variation  of  stress  and  strain  behind  the  unloading 
shock  is  given  by  (5)  and  (7),  showing  that  the  stress 
asymptotically  approaches  zero  and  the  strain  approaches 
the  value  oo^  given  by  (13),  as  shown 

schematically  in  Figure  X. 

The  first  transition  point  occurs  at  Xp  where  the 
elastic-viscoelastic  boundary  first  leaves  the 
characteristic.  The  wave  shape  will  be  similar  in 
form  to  those  shown  in  Figure  X.  However  the  magnitudes 
of  the  initial  elastic  loading  shock,  the  subsequent 
unloading  shock,  and  the  residual  strain 

will  all  decrease.  These  changes  are  shown  schematically 


in  Figure  XI. 


In  the  interval  Xp  to  Xq  in  Figure  III,  the  wave 
form  will  differ  from  those  depicted  in  Figure  X  and  XI 


due  to  the  development  of  Area  H.  In  this  area,  charac¬ 
terized  by  paths  below  in  Figure  IV, 
the  material  is  viscoelastic  with  ^  and  ^ 


decreasing.  This  viscoelastic  unloading  portion  of  the 
pulse  is  bovinded  by  the  elastic-viscoelastic  boundary. 
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The  stress  and  strain  variation  at  a  point  x  between  Xp 

and  x_  is  sketched  in  Figure  XII. 

Ql 

The  second  transition  point  is  x^  where  the  un¬ 
loading  shock  is  absorbed.  As  discussed  earlier,  there 
remains  a  discontinuity  in  the  derivatives  of  CT  E> 

so  that  the  wave  crest  is  cusped  as  shown  in  Figure  XIII. 

It  is  known  from  the  mathematical  theory  of  charac¬ 
teristics  that  discontinuities  in  the  derivatives  of  U 
and  propagate  along  the  characteristic  directions. 

Hence  for  infinitesimal  strain  theory,  the  duration  of 
the  loading  portion  shown  as  Area  E  in  Figure  IV  remains 
constant  and  of  duration  t^,  in  contrast  with  experi¬ 
mentally  observed  pulse  shapes  in  met&ls  and  earth 
materials.  One  way  to  modify  the  theory,  so  that  the 
loading  zone  increases  in  duration, is  to  require  the 
equations  for  the  characteristic  directions  to  contain 
the  dependent  variables.  A  physically  plausible  method 
of  accomplishing  this  is  to  include  the  second  order 
terms  in  the  strain-displacement  law,  thus  making  the 
transition  to  finite-strain  theory. 

8 

As  discussed  in  the  previous  rwftort  ,  the  transition 

* 

to  finite  strain  theory  requires  that  an  unloading  shock 
wave  propagate  along  the  characteristic  emanating 
from  the  point  x  =  0,  t  =  t  ,  with  variable  velocity 
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depending  on  the  strain  in  the  region  into  which  it  is 
noving.  To  analyze  the  variation  of  shock  intensity 
when  the  characteristic  directions  are  variable  and, 
in  general,  not  parallel  to  the  characteristic  direc¬ 
tions  at  the  foot  of  the  shock,  we  assume: 

1.  The  stress-strain  law  during  the  loading  cycle 

-  -Ffck)  -gCCrJ 

2.  The  stress-strain  law  during  the  unloading  cycle 

U"(53  -  O'  -  h 

subject  to  the  condition  (14) 

^  C  CE.ao) 

? 

3.  The  internal  energy  during  the  loading  cycle 

ZZ  /,  7^/3  cx(  (15) 

4o  The  internal  energy  during  the  unloading  cycle 


5o  The  first  law  of  thermodynamics 


cJ  ^  —  ~~  0~  d  oL  0/ 


(17) 
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With  these  definitions  and  assumptions,  the  well- 


known  Hugoniot  function  becomes 


Subcripts  1  and  2  refer,  respectively,  to  points  at  the 
top  and  bottom  of  the  shock.  Write  the  variables  at  the 
foot  of  the  shock  in  terms  of  the  jumps  in  stress  and 
strain  and  thi^ir  values  in  front  of  the  shock: 

oi  ^  o~,-  A(r 

-  Acy 

\/i  -  w  -  A  r 

Then  (18)  becomes 


Write  A  cy  in  terms  of  A  Q~  using  the 

mechanical  shock  conditions  remembering  that  the  shock 
velocity  is  a  function  of  ; 


A  Cx{ 


cr 


(21) 
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Hence  (20)  can  be  written  as  an  implicit  function  of  ^  O' 
and  known  quantities  at  the  top  of  the  shock: 


-  ^ 


<22) 


It  is  interesting  to  note  the  restrictions  placed 
on  the  internal  energy  functions  if  experimental  evidence 
shows  the  unloading  shock  is  absorbed.  If  we  place  , 

(22)  reduces  to 


(23) 

where  ^  and  ^54,  are  not  constants.  This 

requires  that  £  and  be  identical  in  form 

over  some  range  of  (J~  and  oicT  that  must  be 
determined  experimentally.  Over  this  range,  one  func¬ 
tion  will  represent  the  internal  energy  independently 
of  the  equation  of  state  during  different  portions  of 
a  complete  stress  cycle. 

We  illvistrate  the  numerical  calculations  neces¬ 
sary  to  determine  the  jumps  in  stress,  strain,  and 
particle  velocity , using  (15)  and  (16)  to  define  the 
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internal  energy.  Thus  (22)  becomes 


^cr  =c? 


w))ere 


fcYf,) 


/ 


_ 


/ 


The  quantities  A^,  A^j  A^  are  evaluated  in  terms  of 
the  known  stress  and  strain  in  front  of  the  shock.  Since 
(24)  has  two  roots,  that  root  must  be  chosen  which  is 
physically  compatible  with  (19),  the  other  being  extran¬ 
eous  . 

Application  of  (24)  along  the  shock  path  determines 
the  variation  of  stress,  strain,  and  particle  velocity 
along  a  path  in  the  x-t  plane.  The  unloading  problem  is 
completed  by  solving  (14)  simultaneously  with  the  laws 
for  conservation  of  mass  and  momentum,  subject  to  the 
boundary  condition  (j  =  0  at  x  =  0  and  satisfying 
the  previously  obtained  values  of  stress,  strain,  and 


(25) 
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particle  velocity  at  the  foot  of  the  shock  wave.  This 
analysis  applies  only  when  the  shock  path  is  not  a  char¬ 
acteristic  in  the  unloading  region.  The  problem  is  now 
reduced  to  a  new  and  simpler  boundary  value  problem 
with  known  methods  of  solution.  Numerical  solutions 

9 

can  be  obtained  using  procedures  described  by  Lister 
for  one-dimensional  compressible  fluid  flow.  These 
procedures  are  outlined  in  Appendix  I. 
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SUMMARY 


Qualitatively  correct  changes  in  shape  of  a  rec¬ 
tangular  stress  pulse  applied  at  the  surface  of  a  non¬ 
linear  viscoelastic  material  are  obtained  for  the  special 
case  of  infinitesimal  strain  theory  and  a  linear  unload¬ 
ing  law, such  that  the  characteristic  directions  during 
the  entire  stress  cycle  are  unchanged.  It  is  found  that 
the  limitation  of  infinitesimal  strain  does  not  permit 
pulse  broadening  during  the  loading  portion  of  the 
stress  cycle,  which  is  at  variance  with  experimentally 
observed  stress  waves  in  soils.  The  transition  to 
finite  strain,  allowing  for  pulse  broadening,  compli¬ 
cates  the  unloading  theory  so  that  the  mechanical  shock 
conditions  are  not  sufficient  to  determine  the  changes 
in  stress,  strain,  and  particle  velocity  across  an 
unloading  shock.  By  specifying  the  dependence  of  the 
internal  energy  on  stress  and  strain  for  the  material, 
it  is  possible  to  construct  the  Hugoniot  function  which, 
coupled  with  the  mechanical  shock  conditions,  is  suffi¬ 
cient  to  determine  the  transition  across  the  unloading 
shock. 

Once  the  Hugoniot  function  for  the  material  is 
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defin«d,  it  is  possible  to  analyie  the  unloading  shock 
for  quite  general  unloading  laws.  Finally,  the  numeri 
cal  procedures  for  completion  of  the  analysis  are  out¬ 
lined  in  a  form  suitable  for  programming  for  digital 
computers . 
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APPENDIX  I 


We  outline  the  computational  procedures  for  calculating 
the  variation  of  stress,  strain,  and  particle  velocity  behind 
the  unloading  shock  wave  for  the  finite  strain  case.  The 
necessary  modifications  to  reduce  to  infinitesimal  strain 
theory  are  discusmed  at  the  end. 

The  equations  to  be  solved  are; 


(27() 


(28) 


'/x  -  £r  - 


(29) 
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Using,  the  techniques  presented  in  the  previous  report  , 
the  set  of  characteristic  equations  and  directions  corresponding 
to  (26) -(29)  are  found  to  be: 

C4-:  -  J 

dt  '<  J> 

C_ :  ^ 

dt  V  y-> 

ol  cr  V-  p  J_  !/  ^ 

J  f  (33) 

Cot  ^  o 


(30) 


(31) 


(32) 
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crfs)-  (T  = 


(i<) 

i 

Where  S  is  a  parameter  specifying  position  along  the  shock 
path. 

The  boundary  conditions  to  be  satisfied  are: 

(a)  Along  x  =  O,  C7~  =0 

^  '^S=o 

(b)  Along  the  unloading  side  of  the  shock  path 

take  on  the  values  0~(S)^  cries') ^  E.  C  S  ^(5)  which  have 

been  previously  calculated  by  the  continuous  solution  of  (24). 

We  are  specifying  the  values  for  all  the  dependent 
variables  along  this  path.  For  a  solution  to  exist  under 
these  conditions,  we  require  that  the  slope  of  the  unload¬ 
ing  shock  path  in  the  x-t  plane  measured  by  dt  be  every- 

3x 

where  less  in  absolute  value  than  the  slope  of  any  charac¬ 
teristic  for  the  unloading  system  which  crosses  it.  The 
inequality  associated  with  (14)  is  sufficient  to  satisfy 
this  requirement.  Under  these  restrictions,  the  set  of 
characteristics  in  the  unloading  zone  is  shown  schematically 
in  Figure  XIV. 

Equations  (14)  and  (30)-(34)  as  written  are  not  in  a 
suitable  form  for  numerical  calculation.  We  rewrite  them 
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in  incrsmantal  form  after  introducing  the  following  notation 
where  the  relative  position  of  points  P,  Q,  R,  S  are  illustrated 
in  Figure  XVI. 

G  (P)  “  J  ( ‘•‘■(i)-oCiO  •f-L.')  evaluated  at  any  point  P. 

“  (Q)  =  0J(^)  -  oC-c-Cs)  evaluated  at  any  point  Q. 

C  (q)  *  1/2  /5(P)  +  G{Q)J 
D  (g)  =  1/2  /R(P)  +  H(Q^ 

The  last  two  expressions  are  approximations  of  the  average 
value  of  G  and  H,  respectively,  along  any  line  segment  PQ. 

With  this  notation,  the  stress,  strain,  and  particle  velocity 
at  a  point  P  can  be  found  in  terms  of  the  known  stress,  strain, 
and  particle  velocity  at  points  Q,  R,  S  by  the  simultaneous 
solution  of  : 
along 

(X(p)-f  C(0\/(f)  ^  Ofs;  /cs;  (35) 

along  C_ 

'Tfp;  4-  f  dif)  l{p)  ::  (Jld))  -P  V((!))  (36) 
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along 


-h.i.  rjf,- tp).±i.f)= 

The  location  of  P  and  R  are  known  by  prechoosing  the 
grid  made  up  of  increments  ^  ^  and  ^  •  The  location 

of  Q  and  S  must  be  found  as  part  of  the  solution. 

The  calculation  begins  at  a  point  P  at  time  t^  +  A'^ 
on  x=0.  At  P,  0^  ck  j  are  known.  The  point  Q 

is  the  intersection  of  the  C_  characteristic  through  P  with 
the  shock  path.  On'Qe  the  point  of  intersection  is  established, 
all  four  variables  are  known  there  j^rom  previous  calculations. 

Since  the  equations  are  nonlinear,  an  iterative  procedure 
is  used  to  converge  on  the  correct  solution.  In  the  follow¬ 
ing  equations  a  bracketed  superscript  refers  to  the  order 
of  the  iteration.  Nonsuperscripted  variables  do  not  change 
during  the  iteration  procedure  at  a  given  point.  An  iterative 
scheme  for  locating  Q  in  Figure  XV  is: 


(38) 
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i. 

'  «h.r.  C%)-  iSiH) 

I.  For  n  s  1  an  initial  estimate  of  t(^)(Q)  yields  a  value  G(Q) 

(2) 

I  from  which  t'  '(Q)  can  be  calculated.  The  iteration  is  con¬ 
tinued  until  t^  \q)  —  where  6.  is  a  prechosen 

L  accuracy  factor. 

]  Once  t(Q)  has  been  evaluated  to  the  desired  accuracy,  the 

I 

coefficients  in  (36)  can  be  immediately  evaluated  and  the  equation 
I  solves  for  V(  P) . 

The  calculation  at  a  point  P  not  on  the  boundary  is  more 
complex.  We  illustrate  this  more  general  case  in  Figure  XVI 
where  the  new  P  has  been  moved  one  unit  to  the  right  and 

I 

remains  one  isTT  unit  above  the  shock  path.  The  points  Q,  R,  S 

are  respectively  the  intersection  of  the  shock  path  with  C_  and 

C  through  P,  and  the  intersection  of  the  previous  C  characteristic 
o  o 

with  C4-  through  P.  An  equation  from  which  the  time  t(")(S)  can  be 
evaluated  is 


(r)) 

t(s)  -  tcf)  - 

-here  -  C3,(k) 
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To  locate  the  point  Q  we  assune  the  path  of  the  unload¬ 


ing  shock  is  expressible  as  some  numerical  function 
t  =  or  (p 

The  point  Q  is  found  by  solving  the  implicit  function 

t(p]  -  <ffx7o)) 

where 

The  coefficients  in  (35)  and  (36)  can  be  evaluated  by 
assuming  the  initial  estimates 

(fp '  a  + (/?s)j 

atif)  -  I  L  +■  +  •=(‘‘'(v] 

v'h  =  sL  Aq)  f  ^  Asj] 
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/(O 

These  initial  estimates  are  then  used  to  calculate 
by  simultaneous  solution  of  (35)  and  (36)  and  then  solving  (37)  for 
(P)  , since  CT  (O  is  known.  These  estimates  are  then 

used  to  evaluate  more  accurately  the  coefficients  in  (35)  and  (36); 
the  problem  is  then  resolved  as  described  above.  Processes  of  this 

nature  are  stable  and  solutions  can  be  obtained  to  any  desired  _  _ 

degree  of  accuracy. 

When  the  desired  accuracy  has  been  achieved,  the  calculation 
moves  to  a  new  P  point  one  A  )<.  '  unit  to  the  right  and  remaining 
one  unit  above  the,  unloading  shock  path.  The  problem  of 

evaluating  the  stress,  strain,  and  particle  velocity  at  this  new 
point  is  identical  to  the  previously  described  problem  and  the 
entire  process  is  repeated.  In  this  manner,  the  numerical^solution 
is  obtained  along  a  line  parallel  to  the  original  unloading  shock 
path  as : far  out  as  desired.  This  new  path  is  then  used  as  the  base 
onto  which  another  strip  of  solution  can  be  added. 

The  above  procedure  was  outlined  for  the  more  general  finite 
strain  theory.  For  many  purposes,  infinitesimal  theory  will  suffice. 
The  above  process  is  readily  modified  to  include  the  infinitesimal 
strain  case  by  formally  setting' and  its  derivatives  identically 
zero  wherever  they  appear  and  considering  cxi.  as  the  infinitesimal 
Strain.  Under  this  reduction,  (29)  becomes  meaningless  and  should 
be  ignored. 
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APPENDIX  II 


Malvern's  derivations  of  (1)  and  (2)  are  straightforward. 
We  repeat  them  here  since  his  derivations  are  unpublished. 

A  characteristic  equation  valid  along  at  the  top  of 
the  shock  is 

JoT-pC  ,  (41) 

If  the  material  at  the  foot  of  the  shock  is  elastic  with 
Young's  modulus  b/a,  the  characteristic  equation  correr- 
sponding  to  (1)  is 

a.  1/2-  —  K  ”  constant  (42) 

The  mechanical  shock  conditions  yield 

From  (42)  and  (43), 

^  a-  ^  f  A 

^  (44) 

Taking  the  differential  of  (44)  along  C^,  one  obtains 

d  (45) 

Substitution  of  (41)  in  (45)  gives 
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Equation  (1)  is  then  obtained  by  itit^grating  along  from 


t  at  X  s  0  to  t. 
o 


To  obtain  (2),  the  characteristic  equation  at  the  foot 
of  the  shock  is  changed  from  (42)  to 


(46) 


Now  subtract  (46)  from  (41)  to  obtain 

Rewrite  (43)  in  differential  form 

Solve  (48)  for  d  l/^  and  substitute  into  (47)  to  obtain: 


(50) 


Collecting  terms  and  using  cj  OJ one  obtains 


from  which  (2)  may  be  obtained  by  integration  along  from 

t  Tt  o  t . 
o 
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